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ABSTRACT 


The  case  of  a  screamtube  of  cold  gas  gotng  through  a  strong  normal 
shock,  wave  is  treated.  Chemical,  reaction  is  allowed  to  occur  downstream 
of  the  shock,  where  the  pressure  and  the  distribution  of  certain  turbulent 
correlation  functions  are  prescribed  as  a  function  of  distance.  It  is 
shown  that  the  turbulence  modifies  the  flow  field  mainly  through  the  ef~ 
feccs  it  has  on  the  reaction  rates  leading  to  a  hotter  gas  than  if  just 
the  mean  values  of  the  thermodynamic  functions  for  the  turbulent  flow 
were  considered*  The  Gibbs*  tunction  for  the  fluid  is  also  modified  by 
the  turbulent  fluctuations,  thus  leading  to  different  final  equilibrium 
Condi tions . 
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1.  iHTRODUGTimi 

The  preterit  work  was  aoCivated  by  an  attempt  to  find  theoretical 
explanations  of  the  observable  pheftoaena  regarding  trails  left  in  the 
atfflosphere  by  objecl?  flying  through  it  at  hypersonic  speeds. 

-  3 

In  the  case  of  equilibriuia  flew,  it  has  been  shown  by  the  author  (1) 
and  others  (2)  that  for  laminar  flow,  the  trails  left  by  blunt 
objects  could  be  miles  long  depending  on  the  detector  used  by  the 
observer.  It  was  also  recognized  (1),  and  recently  quantitatively  sub* 
stantiated  by  Lees  and  Hromas  (3),  that  when  the  trails  beccsse  turbulent, 
the  cooling  occurs  much  more  rapidly  than  in  the  laminar  diffusion  case. 

It  is  also  known  that  when  certain  metallic  pellets  are  fired  at  high 
speed  or  when  small  meteors  disintegrate  In  the  atmosphere,  they  leave 
luminous  trails  which  scsBetimes  are  longer  than  estimated  frcs  tne  laminar 
or  turbulent  thermodynamic  equilibrium  models  {for  pure  air,  for  example) 
studied  in  the  already  mentioned  references.  This  effect  has  been  easily 
’’explained  away"  in  the  past  by  saying  that  the  ablation  products 
contamiriate  the  high  temperature  air,  thus  leading  to  large  increases  in 

radiation. 

Although  the  foregoing  explanation  may  be  the  correct  one  in  many 
physical  situations  of  interest,  the  question  was  investigated  by  the 
author,  sometime  ago,  regarding  the  possibility  of  the  existence  of 
a  new  effect  caused  by  the  possible  non-linear  interaction  of  the  turbulent 
oscillations  (in  a  compressible  reacting  fluid)  with  the  reaction  rates 
themselves.  Could  this  interaction  affect  the  rate  of  cooling  of  the 
turbulent  wake  of  a  blunt-nosed  body?  The  answer  turned  out  to  be  in  the 
affirmative,  but  the  effect  was  a  function  of  the  magnitude  of  the  oscillating 
quantities  relative  to  their  mean  values.  Some  of  these  quantities  have 
been  measured  in  relatively  low  speed  flows,  and  for  those  cases,  they  are  of 

^Humbers  in  parentheses  denote  references  at  the  end  of  this  paper. 
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small  magnitude.  However^  SletCery  and  Clay  (4)  have  recently  shown 
experimentally  that  the  oscillations  in  hypervelocity  wakes  can  be 
large.  It  was  this  last  result  that  motivated  the  present  quantitative 
study. 

2.  EFFECT  OF  TURBULENCE  ON  THE  CONSERVATION  EQUATlOuS  OF  A  REACTING 

FLUID 

In  order  to  bring  to  the  forefront  the  coupling  of  the  turbulence 
with  the  chemical  kinetics^  we  will  ne‘^lecc  all  normal  gradients  when 
compared  with  the  ones  in  the  direction  of  motion.  Thus^  Reynolds  and 
viscous  Stresses  and  laminar  and  turbulent  conduction  disappear  from  the 
continuity^  momentum  and  energy  equation.  In  general^  this  is  a  very 
important  restriction,  but  is  not  important  for  the  purpose  of  this  paper. 
The  problem  will  thus  be  reduced  to  the  invtscid  adiabatic  flow  in  a 
streamtube,  with  a  prescribed  pressure  as  a  function  of  distance.  However, 
the  conservation  equations,  so  drastically  simplified,  contain  quantities 
which  are  in  general  functions  of  time.  When  mean  values  in  time  are 
taken,  certain  correlation  functions  appear.  The  most  important  effect, 
however,  occurs  in  regard  to  the  modification  of  the  chemical  kinetic 
races . 

2.1  Tne  Conservation  EguattQns 

Tne  conservation  equations  for  a  compressible  reacting  fluid 
with  no  normal  gradients  subject  to  a  prescribed  pressure  gradient  will 
be  written.  The  reason  for  treating  essentially  the  case  of  a  streamtube 
which  is  Subject  to  a  given  pressure  distribution,  is  that  it  is  the 
simplest  case  that  could  be  considered  in  which  the  main  features  of  a 
hypersonic  wake  of  a  blunt  body  can  be  kept  in  the  problem^  the  strong 
pressure  gradients  which  govern  the  chemical  kinetics  of  the  problem. 

This  streamtube  could  be  interpreted  as  the  first  approximation  to  the 
center  streamtube  of  a  turbulent  reacting  wake.  In  order  to  keep  the 
problem  simple,  the  analysis  will  be  restricted  to  a  diatomic  gas  which 
can  dissociate  and  recombine.  I.et  any  scalar  quantity  time  varying 
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quanciUy  q  be  decopsposed  into  q+q*,  where  q  represents  the  tise  averaged 
value  and  q'  is  the  time  dependent  parr.  Also,  let  bars  over  a  quaritity 
denote  a  time  average  taken  at  a  given  point  in  space  and  will  be 
used  only  when  its  deletion  may  lead  to  ambiguity.  The  equations  for 
the  mean  quantities  in  one'diffiensionai  flow  are  given  below  ^5) 

/  \ 

H 


Momentum: 


du 

Pu  -7“ 

I  dx 


+  0‘u‘ 


du 

dx 


Q 

dx 


where  u,  p,  p  and  T  denote,  respectively,  the  time  averaged  values  of 
velocity,  pressure,  density  and  temperature,  u'^  is  the  Reynolds  stress 
in  the  x-dlrection,  s  denotes  axial  distance,  non-dimensionalized  for 
example  with  respect  to  body  radius  r^;  l.e.  s  «  ^  undissociated 

{cold)  gas  constant,  the  primes  denote  time  varying  quantities,  and  the 
bars  denote  time  averages  at  a  fixed  position  in  space.  Ute  Subscript 
CO  denotes  values  of  the  quantities  in  the  undisturbed  gas  ahead  of  the 
projectile j  Our  streamtube  will  consist,  at  the  extreme  left,  of  the 
free*8tream  gas  which  first  goes  through  a  normal  shock  wave  (discontinuous) 
across  which  no  chemical  reaction  takes  place.  Downstream  of  the  shock 
wave,  chemical  reaction  is  allowed  to  proceed  with  a  prescribed  ressure 
niscory  which  can  be  approximately  prescribed  a  priori  for  a  particular 
body  shape  because,  for  our  purpose,  it  is  not  a  sufficiently  sensitive 
function  of  the  reaction  rates. 


Continuity: 


where  A  denotes  the  cross  section  of  the  streamtube,  and  the  subscript 
CD  denotes  free  stream  conditions. 
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where  h  denotes  enthalpy. 


From  statistical  mechanics,  if  it  is  assumed  that  the  vibra¬ 
tional  energy  is  always  in  equilibrium  with  the  translational  energy 
(i.e,  vibrational  excitation  rates  are  infinitely  fast),  and  this  is 
not  a  significant  restriction  for  the  present  purpose,  the  enthalpy  of 
a  dissociated  diatomic  gas  cai  be  written  as  a  function  of  temperature 
and  atomic  mass  fraction  y,  as 


RT 


CO  00 


(5  +  ^  )y  +  <  I  +  p)  (1-y) 


1  e. 


J  P 


-  (1-y)  p  Qxp  (e^/T^  jl  +  Y  ^  ^  2  ^  ‘  ^ 

I  1 

where  9  is  the  vibrational  constant  of  the  molecule  and  ot  m  5^365°K 


rr 


9  /T 

V 


exp  (e,/T)  -1 


. . . . . . . 


where  and  stand,  respectively,  for  the  dissociation  and  recombina¬ 
tion  terms  in  the  kinetic  equation,  and  M  is  the  atomic  weight  of  the 
chemical  species.  The  novel  results  that  will  be  presented,  will  depend 
heavily  on  Eq.  6  which  -•ill  be  derived  in  detail  in  the  next  section. 

The  unknowns  are  A,  u,  y  and  T  as  functions  of  s,  as  well  as 
all  the  correlation  functions.  l^\ese  will  have  to  be  specified  as  a 
function  of  s  and  will  be  discussed  in  Section  3. 

2.2  The  Chemical  Kinetic  Equations 

There  is  no  particular  virtue  in  writing  v»hat  follows  in  general 
form  for  any  diatomic  molecule.  If  done  in  general,  some  of  the  expressions 
will  become  too  lengthy  and  cumbersome.  Any  reader  ir.erested  in  using 
the  present  approach  for  his  particular  reaction  will  certainly  want  to 
rederive  all  the  expressions  on  his  own.  Thus,  we  will  specialize  all 
the  following  discussion  to  oxygen  as  a  typical  component  of  air. 

Consider  the  reaction 


0  +  U  +  M 


0^  +  M 


where  the  third  body  M,  could  either  be  0  or  O^.  The  k’s  are  the 
reaction  rate  constants.  Let  a  superscript  on  them  indicate  the  type 
of  third  body  under  consideration,  and  also  let  |0  J  denote  the  instantaneous 
concentration  of  oxygen  atoms  in  moles/cm  ,  the  reaction  rate  can  be 
expressed  as 


-  [oj  -  <  w  -  y  ra 


The  concentrations  can  be  written  in  terms  of  mass  friction.,  as 


moles 


(g3»/mole) 


RR-5,  3-62 


where  gs/sale  and  p  is  the  density  of  the  gas  aixcure. 

Froffi  &ef.  (6) 


.  Oo  A.  0 

■"d  ^  ^  -  25k,  , 


where 


k . 


3.6  X  10‘«  T-1  expC-  ^ - 

T  laole  sec  * 


>  10 


and  T  is  in  K,  We  also  know,  fros  theoretical  chemical  kinetics, 
'^iihin  the  restriction  that  a  Soitzciann  distribution  exists  at  the  non* 
equilibriiffl  tesperature,  that 


,  ,  0*1  1 0 

k .  k,i  k, 

k  ,  0*  ,0  ^ 

r  k^ 


0-0 


0<jikli.xbrxuni  constsnc*  Fox’  t6sip€rdCyX‘€5 
O'^c^T  in  oxygfin  w?isn  flying  up  vo  veiociui^s  of  sl>outi  23^000 


given  frcfe  statistical 


S„CT)  -  2.93  X  10*^  (1 


SccnariitS  DV 


^  r\ 


-228 

H-  3  e  ^  -f-  e 


J25.9  \T1 


I  3  +  2.027  e 


-U^ 

1 


-59.365 
e  T 


soies 


In  tne  range  between  300C  K  and  8ijuO*^K,  the  above  expression  can  be 
approxisated,  witn  less  than  lO^o  error  (6}  by 


Kit) 

c 


1.2  X  10^  expC-  — y. 

cm 
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Thus#  ffos  Eqs,  10#  ll# 


and  13 


k®2  -  9k^ 


kj  -  25k^ 
6 


and 


k  ■  3  X  10^^  — T” — 

r  issole*  sec 


Inserting  Eqs,  9#  10  ana  14  into  S  yields 


4^  «■  k,  P  (20.5  y  4-  4.5)  -  k 
lie  ®  I  * 


'4iy  +  9)# 


where  the  tilde  has  been  added  to  eisphasize  that  the  quantities  are#  in 
the  case  of  turbulent  flow#  coesposea  <',,3  tiise  averaged  and  an 
oscillating  value. 

g^writiag  Eq.  15  in  terss  of  sean  tnO  cisse  dependent  quantities  gives 

^  i  (k^  -i-  k*)(pi-K*)(l-y-y’)  1^0.5  (y+y')  4.^ 

dt  dtW  dCil  I—  — * 


0 

M" 


where  we  have  used  the  fact  that 


(y+y') 


I4i(y47’) 


y'(t) 


and  y(t) 


yCt)  -  y<t) 


y(c). 


Carrying  out  the  operations  in  Eq.  16#  taking  tise  averages#  and 
considering  that  the  tise  average  of  a  prised  quantity  or  Its  derivative 
is  zero#  leads  to 


dv 

dt 


where 
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Iji“(^*5+20.5y)  (l-y)  (pk^+l^^kp  +  py*k^  (16-4Iy)  k^p'y*  (l6-4ly) 


■20.5  pk^  y 


+p'y*k*.  (16-4 ly) 


20.5Cpy*‘k! 


+  k 


7^) 


•20.5  y'^p’k* 

r  ° 


p  _  _ 

\  i£(Q+9)k^py  p*¥*  -r  2P  py'p'k^  +  (Q+9>p^y  y 'k^4Pk^y‘p •24<^k^*y  ^ 


’»^  *  2+ok  £>^T^ 


+  (Qf9)k^y  P’^y*  -i-  2Qk  C^'^p.  ^  2  (Q+S)  py  P*  >' ’  k  ’  +  Py^  k*p*^  +  Qf^kjy^<f4lk 
*  *  ri  I  ;  ;  r  r*  i  r  ri 

2  - ? 

'*^rf  ^  *r  P  ‘‘r  ^ k^y'p' 
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♦  [q. 


t3i-kt  _t_  AT ...2  ,3, 


Qk  p*  y’‘  -I-  2Fy0*2y«k*  +  2Qpy*"p‘k‘  4- 41y^  k'y’p'‘  +  82k  py‘^p‘  +  4ip‘y‘"k 
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X-3 


^y  +  9. 


Q  «  123y  +  9. 


-ne  first  ters  in  £qs,  18  and  ly  is  the  usual  one  that  appe-^*^  ^ 


cars  t»hen  the 


iiuw  IS  iasinar,  without  r.avin 


ing  sny  experimecitSi  information  on  the 

u^H^vior  of  tno  corroi^tion  iorivt^ons  thst  sppe^r  m  tnese 

It  could  De  conceivable  that  the  ters  I  which  in  laminar  flow  is  the 

c 

contribution  to  the  net  rate  due  to  dissociation  say  here  become  a 
recosbination  term,  A  similar  cosusent  could  be  made  about  Mhen 

such  a  thing  nappens^  the  phenomenological  approach  that  led  to  stating 
the  cheslcal  kinetics  in  the  form  of  Eq.  8  may  require  modification, 
aquation  18  contains  up  to  fourth  order  correlation#  while  Sq,  19  contains 
up  to  the  sixth  order.  It  should  be  noted  that  in  deriving  Eqs*  18  and 


l9  no  terss  vore  neglected^  and  that  no  restrictiva  assifflption  was  issde 
regarding  the  sire  of  primed  quantities  relative  to  the  corresponding 

®ean  values  (i*e,  p*/p)- 

I  i 

In  order  to  reduce  Eqs.  18  and  19  to  a  tractable  form,  it  will 
be  necessary  to  derive  expressions  or  evaluate  the  correlation  functions 
that  appear.  Very  few  correlations  higher  than  the  second  have  been 
made  in  turbulent  fluid  mechanics  which  are  applicable  to  the  present 


problem.  If  the  oscillating  quantities  are  small  compared  to  their 
respective  mean  values,  correlations  higher  than  the  second  will  not  be 
of  importance.  In  the  case  of  hypervelocity  wakes  Slattery  and  Clay  (4) 
show  that  0.5,  which  is  far  from  small.  However,  in  order  to 

make  the  problem  tractable,  ■it  will  be  assumed  in  the  present  paper  that 
only  second  order  correlations  are  of  importance.  After  obtaining 
numerical  results  based  on  this  assumption,  they  will  be  examined  and 
evaluated  in  terms  of  it.  From  £q$.  18  and  19  the  second  order  correla* 
tions  required  are 


I 


y  -4. 


^  *1.  • 
r  r- 


and  » 


2  ,2 

These  functions  should  be  expressed  in  terms  of  p’  ,  p’T*,  and  T*  ,  It 

I  f  ^ 

will  then  be  necessary  to  expand  each  of  the  prime  quantities  in  Eq.  21 
in  p*  and  T‘,  take  products  and  then  time  averages. 

From  the  perfect  gas  law 

p+p'  -  R(p+P’)CT^T‘)(l-HH-y*) 
or 

p+p‘-  R[pT+(|^)]]  a+y-^y‘)  [220 


Since  the  perturbation  pressure  will  depend  on  the  square  the  perturba¬ 
tion  velocity,  it  will  be  assumed,  as  a  first  approximation,  that  p'*  0. 

A  discussion  of  this  point  is  given  fay  Kistler,  with  regard  to  compressible 
boundary  layer  work,  in  Ref.  7,  p.  295.  Simultaneous  measurements  of  all 
the  oscillating  quantities  in  Eq.  22  have  not  been  made,  ev^en  when  y+y*  *  0. 


Sose  pressure  seasurements  made  in  boundary  layers  by  different  investi¬ 
gations  for  various  flow  conditions  are  briefly  reported  by  Willssarth 
in  (Ref*  8,  p*ll2)  from  which 


"11/2 


10 


P  2 

where  y  is  the  ratio  of  specific  heats,  and  M  is  the  Mach  number  of 
the  flow  external  to  the  boundary  aajfer.  In  order  to  apply  this 
result  to  a  turbulent  wake,  it  is  thought  that  the  significant  Mach 
number  should  be  based  on  the  velocity  defect  across  the  turbulent  cure 


of  the  wake-  for  a  blunt  object  Eq.  22  will  approximately  be 

reduced  to  !  p  *^1  *  10 


Setting  p*«»0  in  Sq*  22a  and  solving  for  y‘  gives 
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or^  if  using  Eq.  22b 
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The  time  averaged  terms  that  appear  in  Eqs.  24a  and  b  are  of  higher 
order  than  the  others,  and  can  therefore  be  neglected.  In  order  to  end 
up  witn  oscillating  quantities  of  a  single  variable  in  each  term,  and 
since  the  oscillating  products  are  such  that  are  not  any  easier  to 
measure  than  the  single  quantities,  we  will  expand  the  last  term  in 
parenthesis  and  the.  denominator  oi  Eq.  24a  into  a  power  series  of 
p*/p  shd  T'/T,  l»hen  the  dame  is  done  for  kj  and  k*  the  results  are 

I  i  dr 
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,  and 


Taking  the  products  necessary  to  construct  the  expressions  in  Eq.  21,  and 
averaging  in  time  yields  “ 

«  I  >n  i  m  f  fcN 

riple  correlations 


/p'2  .  . 

pij'  T'"^^ 

I  i  2  ^ 

p*y’  ^  -(1+y) 


triple  correlations 


oT 


q*K  •  A 


^  +  triple  correlations 

d  I  a  D* 
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D  i  ,  T 
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Inserting  Eqs*  26  into  Eqs.  18  and  19  yields^  after  dropping  third  and 
higher  order  correlations^ 


I 


d 


too  P 

ICO 


I  2 

f 


+  f. 


where 


(4.5  +  20.5y)(l-y), 

(l+y)  |^(16-41y)  +  20.5(l+y)j  , 

(l-Ky)(36.5  20. 5y), 

(4.5  +  20.5y)(l-y)  q-(l+y)  {  q(16-4ly)  +5^, 


and 


CO 


1 — I  k 


T*  <  a  •t*  • 

U  +  g.  —  +  g^  ^  ^  ^  - 

T  P 


h  bt 


{ 


wnere 


(41y  +  9)y\ 


g,  •  3 (l+y) 


y 


“  (4ly+6)y  +  (41y+3)  (1-K') 


.T1 


(4iy+9)y+3(l+y)  (41y4-3)  ^^^^-2(4Iy-6) 


r. 


8^  *”2  C41y+6)  (l-fy)y+2  (41y+3)  (1+y)^  -  ~  (4ly+9)y^  . 


It  should  be  pointed  out  that  and  g^^  in  Eqs.  28  and  30^  are  the 
usual  terms  that  appear  in  the  chemical  kinetic  equations  when  there  is 
no  turbulence.  The  remaining  f's  and  g’s  can  be  interpreted  as  influence 
coefficients  of  the  correlation  functions  on  the  reaction  rates, 

Me  will  discuss  next  how  to  arrive  at  some  reasonable  values  for 
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the  correlation  functions  necessary  to  solve  the  probles. 

3.  ESTIMATES  OF  TOl  GORRELATIOH  TiJliQTlCMS  BASED  m  S(m.  EXPEIlMEl«!AL 

BATA 

There  are  many  different  correlation  functions  that  appeared  in  the 
foregoing  discussion.  Only  few  of  thes  have  ever  been  measured  regard¬ 
less  of  the  fluid  mechanical  problem  under  diacussion*  Therefore,  based 
on  whatever  experisaental  data  are  available,  we  will  have  to  at  least 
evaluate  the  order  of  magnitude  of  the  functions  involved.  In  order  to 
do  this  in  a  reasonable  manner,  we  will  quickly  and  superficially  review 
in  what  follows  some  of  the  typical  values  of  certain  correlation  functions 
that  have  been  measured  in  different  flow  problems,  i,e,  boundary  layers, 
jets  and  wakes, 

3.1  Measurements  in  the  Boundary  Laver 

Typical  measurements  are  those  presented  by  Klebanoff(9)  of  the 


u  —  10  in  Incompressible  boundary  layers;  those  of  Sis tier (7)  of 


♦  € 

f^jl/2/^  ^1Q-I^fu.2)l/^^^  -7  X  10*^,  (p>^u«0.l4  in 

compressible  boundary  layers  between  mach  numbers  of  1.72  and  4,67  where 
the  stagnation  temperature  was  300^K. 


3.2  Measurements  in  Turbulent  Jets 

Data  for  incompressible  jets  can  be  found  in  Townsend's  book  (10). 
The  maximum  disturbances  occur  near  the  center  of  the  jet  and  gives 

-0.15,  u'v'/u^~  .008.  Corrsin  and  Uberoi  (11)  measured  for 
a  heated  jet  T  —  10 

3.3  Measurements  in  yakes 
3.3,1  Incompressible  Wakes 

For  incompressible  wakes  flows  Townsend  shows  that 


-  .006, 
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3.3.2  Compressible  Makes 

Demecrlsdes  (12)  sseasured  (pu)  ’/p^  *0.40  in  che  wake  of  a 
cylinder  at  free  stream  mach  number  of  5.8  and  Reynolds  numbers  between 
18^500  and  48^500.  Slattery  and  Clay  (4)  have  taken  schlieren  photographs 
of  the  wake  left  behind  a  spherical  projectila  of  1/2  inch  in  diameter 
flying  at  8370  ft/sec  in  air  at  a  pressure  of  30  mra  of  mercury,  'Riey 
took  densitometer  tracings  of  the  negative  at  several  portions  downstream 
of  the  projectile  and  derived  from  them  the  density  fluctuations  as  given 
in  the  following  table: 


x/d 

lu 

aO 

400 

800 

0.05 

1— _ 

0.25 

0.90 

0.50 

The  fluctuations  in  density  can  be  seen  to  be  quite  large^  and  they  are  the 
first  experimentally  derived  values  obtained  at  hypervelocities.  These 

values  should  be  compared  with  p’u‘  /pu  *0,4  obtained  by  Demetriades  (11). 

I  I 

Most  of  the  oscillation  of  the  mass  flow  correlation  must  be  due  to  the 
oscillation  in  density,  because  if  it  were  due  to  the  oscillation  in 
velocity,  Che  kinetic  energy  invested  in  che  tv.rbulence  would  be  pro¬ 
hibitively  large.  Before  use  is  made  of  a  high  value  for  the  density 
fluctuation,  an  attempt  will  next  be  made  of  explaining  why  such  a  value, 
although  high,  could  be  justified  in  che  case  of  hypersonic  wakes. 


slder  the  wake  of  a  hypervelocity  sphere  (Fig.  1)  that 
has  a  turbulent  wake,  Lees  ana  Hromas  (3)  have  analyzed, in  some  detail, 
the  behavior  of  such  a  turbulent  wake  in  equilibrium  flow.  Within 
the  turbulent  core,  one  would  expect  to  have  eddies  of  various  tempera¬ 
tures.  Tne  range  of  temperature  variation  of  the  eddies  would  be 
between  the  low  temperature  in  the  invisctd  flow  at  the  turbulent  edge 
and  the  highest  value  in  the  turbulent  core.  Hromas  (13)  has  given  che 
author  some  of  his  calculations  for  a  turbulent  trail  of  a  sphere  of  1/2 
inch  in  diameter,  flying  at  sea  level  at  9500  ft/sec.  At  this  velocity 
there  are  no  chemical  reactions  and  che  enthalpy  on  the  center  line  h  ,  of  the 

O’ 

trail  divided  by  che  edge  value  h-,  is  equal  to  the  corresponding 

temperature  ratio,  T  /T.. 

o  £ 
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For  a  blunt  body,  close  to  t he  wake's  neck,  T^/Tj  Should  be  near  unity; 
very  far  downstr 'a®,  when  all  the  gas  has  cooled,  —>1.  In 

between,  the  outer  Inviscid  wake  cools  rapidly  and  this  ratio  should 
have  a  peak.  fig.  2  shows  this  ratio  for  two  bodies  under  different 
flight  conditions. 

We  will  use  the  9500  ft/sec  curve  to  arrive  at  some  conclusions 
regarding  the  turbulence  structure  necessary  to  make  possible  a  value 

j  2  \  T  /'J 

oi\p*  j  0,5,  as  obtained  by  Slattery  and  Clay  (4).  Let  us  assume 

that  in  the  turbulent  core  we  have  a  mixture  composed  of  lumps  of  hot 
and  cold  gas,  which  if  no  radial  gradients  of  pressure  allowed,  represent 
spots  of  low  and  high  density  gas.  Let  p^,p|^  and  p  represent  respec¬ 
tively  the  density  of  the  cold  gas,  the  hot  gas  and  the  mean  value* 

We  will  assume  that  the  structure  is  periodic  in  space  with  a  period 
of  length  L.  We  will  neglect  that  the  gas  may  be  reacting  chemically. 
Thus 


The  magnitude  of  the  temperature  for  the  cold  gas  would  not  be  lower 

(although  it  say  be  higher)  than  T^, the  temperature  at  the  edge  of  the 

turbulent  wake.  T.  would  be  higher  than  the  mean  temperature  calculated, 

n. 

for  example,  by  Lees  and  Hromas,  Thus,  would  be  approximately  given 


Let  the  value  of  the  square  of  the  density  fluctuation  be  given  by 


where  the  total  density  0  is  given  by  p  ■  0+p‘,  and  the  density 

^  'll 

structure  is  assumed  to  be  as  shown  in  Fig.  3.  For  that  structure,  tne  mean 
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Clay's  measurements), is  there  any  hope  to  find  a  value  of  €,  from  Eqs.  34 
and  35,  that  would  make  physical  sense  for  the  model  of  Fig.  3*  Physical 
sense  means  that  the  value  of  €  satisfies 

0  i  00 


should  be  small  compared  to  unity;  that  is,  only  less 


and  even  more,  c 
than  half  of  the  mixture  should  he  cold  gas. 


If  the  value 


of  P, /P  from  Eq.  34  is  used  i..  Eq.  35,  solving  tor  €  leads 
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where  the  §ub?icript5  1  and  2  -denote  the  two  possible  solutions  of  th 
probleHi,  Both  v'alues  of  e  and  p  .  narti-sular Iv  the  second  ones,  are 
reasonable  because  sucn  vaiues  snoum  oe  pnvsxcaily  reaitsabie* 
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Differentiatioa  ot  Eq.  51  with  respect  to  insercion  of  Eq«.  42  and 
50  lead  £c 


la  order  to  solve  the  differential  equations  42  and  52  it  is  necessary 
to  prescribe  the  pressure  on  the  streastube  as  a  function  of  s.  Since 
^»e  are  interested  in  a  blunt  body,  we  will  prescribe  a  Newtonian  pressure 


distribution  on  the 


*05^ 


0^  s^  0.96 (-55°)  ^3J 


wner e  Si^ne  auDauript  s  den-^te^i  stagnation 
pressure  at  90''  away  frosi  the  stagnation 


conditions.  Let  p  ,  denote 

sh 


point,  then  the  pressure  can 


be  expressed  (14)  as 


wnere  "  .  can  be  decersined  (15)  from 
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where  is  Che  free  stream  Hach  number.  Between  55°  and  90°  the 
pressure  can  be  fitted  with  an  expression  of  the  form 
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Finally,  £he  corr elation  functions  accessary  in  Eqs.  11  and  29  are 
pT^/p^  gjj4  p*T*/pT.  The  calculations  will  be  made  for  the 
case  idiere  the  streastube  has  no  turbulence  £r<®  s-0  to  s»3,  it  buiids 
up  linearly  with  distance  until  it  reaches  s»4  (i.e.  in  one  radius)  to  a 
value  that  stays  constant  with  s  (Fig.  4),  i.e. 
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5.  DXSCdSSXC^t  OF  HSHBRlCAh  RESUUS . 

Before  deciding  iAat  values  of  a^,  a^  and  a^  should  have  rn 

the  calculations,  it  is  important  to  realiac  sane  of  the  restrictions 
that  have  to  be  imposed. 

Only  a^  has  been  measured  (4),  Thus  the  other  a‘s  have  to  be  con* 
jectured.  a^^  and  a^  are  positive  definite.  In  a  non-reacting  gas, 
y  -  constant  and  if  in  Eq.  22  a  p  »  0,  p‘t‘  would  always  be  a  negative 
quantity,  and  further  more  would  be  aero.  However,  in  our  calculations 
p ’f *  was  assumed  to  be  negative  or  aero.  Ho  calculations  were  fsade  for 
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p’T'  positive.  It  also  foliovs. 
Various  other  relative  values  of 


troa  the  aon-reacting  ease, 
and  a«  were  used  also. 


that  a^.  .j. 


A  physically  iapossibie  combination  of  a's  manifests  itself  in  an 
obvious  way  by  becoming  negative,  or  buth  have  to  be  positive 
at  all  times, 

«* 

If  the  a’s  are  sept  smaller  than  10  the  present  problem  as  treated 


herein  Is  approximately  correct.  For  larger  values  of  the  a’s  the  triple 

and  hi^er  order  correlation  terms,  which  have  been  left  in  the  calculation, 

become  Important,  These  could  be  kept,  but  the  difficulty  is  to  prescribe 

reasonable  values  to  them  in  a  particular  calculation,  since  not  only  have 

they  not  been  measured,  but  it  may  be  impossible  to  do  so. 

Fig.  4  shows  the  t&sperature  X  and  atomic  oxygen  mass  fraction  y, 

along  a  streamtube  as  a  function  of  downstream  distance  measured  fr^  the 

normal  shock  non-dimensionalired  with  respect  to  nose  radius.  Two  spheres 

have  been  used,  one  of  1/2  inch  T.,4ius  and  the  other  of  5  inches,  flying 

•2 

at  20,000  ft/sec  in  an  oxygen  atmosphere  wtiere  density  is  10  the  normal 
value.  When  the  temperature  reaches  a  high  value  while  y*0. 

After  dissociation  the  teaiperature  drops  and  y  builds  up.  Fig.  4  shows 
tne  History  or  the  rlow  tor  s>  lu  ,  It  can  be  seen  that  for  large  values 
ot  s,  after  the  pressure  has  long  ago  reached  the  free-stream  value,  the 
variables  reach  a  constant  equilibrium  value.  The  reason  for  the  equi¬ 
librium  state  for  each  sphere  Is  not  the  same  is  due  to  different  en¬ 


tropies  of  the  flow  in  each  case:  for  the  smaller  s^ere  the  flow  is 
further  out  of  equilibrium  during  the  expansion  process,  thus  leading 
CO  a  larger  entropy  in  the  final  state,  which  represents  a  higher  final 


temperature. 

It  should  be  noted  that  the  final  equilibrium  state  will  be  affectec 
by  the  presence  of  turbulence.  This  can  be  clearly  seen  from  the  reaction 
rate  equations  when  they  go  to  equilibrium.  One  other  way  of  looking  at 
this  fact  is  that  the  turbulence  terms  modify  Gibbs'  function,  vdtlch  when 
minimized,  leads  to  a  different  equilibrium  condition. 

Figs.  5  through  7  present  some  of  the  results  obtained  for  the  case 
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DISTANCE  DOWNSTREAM  OF  NORMAL  SHOCK  WAVE,  S^x/r^ 
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FIC.  4  T^IRATURE  AHD  AT(M  MASS  FRACTION  DiStRIBU 
TI<^’  OS  AXIAL  STREAMTUBE.  LmiMB,  FLOW  FCR 
two  SPHERES . 


FIG.  5  TEMFERATLRE  DISTRIBUTION  ON  AXIAL  STREAMTLRE. 
A  CQ4FARISa^  BETWEEN  LAMINAR  AND  TURBUmiT 
FL^X^.  SFtffiRE  OF  5  IN®  RADIUS* 
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j-  *  5  Indies.  Xs  these  tigutes  the  curves  should  be  cospered  with  the 
fiQfj=turbulent  results.  Xhe  soire  curves  ere  to  be  trusted^  wntte  the 
clashed  ones  csnnet  be  proven  to  be  correct .  However^  the  feet  thet  they 
form  s  reasonable  family  is  seaningtul.  In  Fig-  5  for  example  it  shoula 
be  noted  that  although  the  final  equilinriua  t^perature  is  increasea  by 
at  most  300®K  dten  the  flow  is  turbulent,  between  s  »  10*  the  Increase 
ean  be  near  lOOO^K*  This  hotter  vase  eouid  then  lead  to  signifteantLy 
more  radiation.  Using  the  data  of  these  tigures^  it  is  possible  to  caleu* 
late  the  absolute  magnitude  of  the  radiation  for  the  different  asstffiee 
situations.  Such  results  would  be  more  sensitive  to  the  turbulence  than 
the  temperature  variations. 

Figs.  8  through  10  give  similar  results  for  the  case  where  the 
sphere  nose  radius  Is  1/2  inch.  It  can  be  shown  that  if  the  correla¬ 
tions  neglected  in  £gs.  42  to  49  were  included  in  the  calculations^ 
higher  temperatures  would  ensue. 

6.  C®ICUIBIKG  REMARKS 

From  the  foregoing  results  it  ean  be  said  that  turbulent  chemical 
reacting  flows  under  the  conditions  studied  can  lead  to  flow  fields^  which 
at  any  given  station  are  of  a  higher  temperature  than  without  turbulence. 
The  fact  that,  in  hyparvelocity  wakes,  the  mean  values  of  the  turbulent 
correlations  ean  be  large,  leads  to  the  possibility  of  significant  effects, 
although  we  have  not  proven  this  with  all  the  necessary  rigor.  The  sur¬ 
face  has  just  been  scratched  and  much  work  remains  to  be  done. 
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